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The possible existence of the Anderson transition in 2D systems with- 
out interaction and spin-orbit effects (such as the usual Anderson model) 
becomes recently a subject of controversy in the literature [U EJ [3]. Com- 
parative analysis of approaches based on generalized Lyapunov exponents 
is given, in order to resolve controversy. 



1. Introduction 

Possible existence of the Anderson transition in the 2D case becomes recently a sub- 
ject of controversy in the literature. Kuzovkov et al [T] studying a growth of the second 
moments for a particular solution of the quasi-lD Schroedinger equation and interpret- 
ing results in terms of the signal theory came to conclusion that the first order Anderson 
transition exists in the usual 2D Anderson model. This incredible result was analyzed by 
Markos et al [2] and the opposite conclusion was drawn: behavior of the second moments is 
qualitatively different from found in pQ and contains no evidence of the metallic phase. On 
the other hand, the present author [3] used the analogous approach and came to conclusion 
on existence of the Kosterlitz-Thouless type transition, which however can be realized not 
in all 2D systems. 

The present paper has an aim to analyze the existing controversy and is organized in 
the following manner. In Sec. 2 we give brief exposition of paper [3], which we consider as 
the proper treatment of the problem. Approach by Kuzovkov et al pQ is discussed in Sec. 3: 
the results for the second moments are shown to be partially correct but their interpreta- 
tion is not satisfactory. Section 4 deals with the approach advanced by Markos et al [2]: 
their statements are shown to be related with improper calculation of the matrix product. 
Conclusions are made in Sec. 5. 



2. Second moments of the Cauchy solution 
and their relation with the Anderson transition 

Consider the 2D Anderson model described by the discrete Schroedinger equation 

+ 1pn,m-l + V n>m ^ n>m = El[> 

n,m V / 

and interpret it as a recurrence relation in the variable n, which we accept as a longitudinal 
coordinate. Initial conditions are assumed to be fixed on the left end of the system, while 
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the periodic boundary conditions are accepted in the transverse direction, ip n ,m+L = VVi,™- 
Site energies V n ,m are considered as uncorrected random quantities with the first two 
moments 

( ^n,m ) , ( Vn,mVn/ ,m' ) ^n,n'^m,m' • (2) 

The growth of the second moments for this problem can be studied using the old idea by 
Thouless |1] based on the observation that variables ip n ,m are statistically independent of 
Vn )T a with the same n. The main quantity of interest is (^J; solving (1) for i/j n +i im and 
averaging its square, we can relate it with the pair correlators containing lower values of 
n. Deriving analogous equations for the pair correlators, we end with the closed system of 
difference equations for the quantities 

•E-m,m' (^0 = ( 1pn,m1pn,m' ) i 

ym,m' (jl) = {lpn,m1pn—l,m' ) ; (3) 
Zm,m'\P') = ( 1pn—l,m1pn,m' ) > 

which for E = has a form [3] 

x m>m '(n+l) = W /2 5 m>m 'X mim /(n)+x m+ i >m / + i(n)+x m _i im / + i(n)+x m+ i im /_i(n)+x m _i im /_i(n) + 
+x m , m '(n - 1) + y m +i,m'(n) + y m -x, m <{n) + z TO , m '+i(n) + ^ m ,m'-i(«) (4) 

ym,m' %m-j-l,m' (p*) •^m— l,m' yl) Zm,m' yl) 
Zm,m' 1) •£m,m'+l(^) •^m,m'— l(^) 2/m,m' v • 

This is a set of the linear difference equations with independent of n coefficients and its 
solution is exponential in n [5] 

The formal change of variable is useful 

■Em,m' = ■Em,m' —m ^ra,l j etc., (6) 

where I — m! — m. Substitution of (6) and (5) to (4) gives the difference equations whose 
coefficients contain no m dependence and their solution is exponential in m 

x m ^ = x l e ipm , etc., (7) 

where allowed values for p are determined by the periodical boundary conditions in the 
transverse direction. Excluding y m j and z m> i from the first equation in (40), we end with 
the equation 

xi +2 e~ ip + xi- 2 e tp + V5i fi xi = exi , x i+L = xi , (8) 

W 2 sinh/3 



e = 2 cosh /3 , V 



cosh (3 — cosp 
2 



describing a single impurity in a periodic chain. The positive exponents (3 S for finite odd 
L are determined by equation^] 



2(cosh/? s - cosp s ) = W 2 coth((3 s L/2) , 



Ps 



2tis/L, s = 0,l,...,L-l. (9) 



Allowed values of p s and j3 s become dense in the large L limit, and the quantities j3 and p 
can be considered as continuous; the minimal value of (3 is realized for p = it and can be 
easily found in the large L limit 

cosh" 1 (W 2 /2 - 1) , W 2 > 4 



A, 



- tanh" 1 (W 2 /A) , iy 2 < 4 



(10) 



21nL - 21nlnL + . . . 
L 



W 2 



4. 



The character of solution is qualitatively changed at W c = 2; this value of disorder corre- 
sponds to a singular point found in pQ but its physical sense is not evident.il 

Interpretation of results can be given using the finite-size scaling approach in the form 
given by Pichard and Sarma [6][1. A quasi-lD system is always localized and the localization 
length £i£> can be introduced for it. The knowledge of this length can be used for the study 
of the Anderson transition in the higher- dimensional systems. It is easy to show that 
£id const in the localized phase and £,id/L — > oo in the metallic state in the large L 
limit [S], E] . If we introduce a scaling parameter 



9(L) 



L 



[ID 



then it increases with L in the metallic phase and decreases in the insulator phase. In the 
framework of one-parameter scaling hypothesis a following relation can be postulated [7] 



g{L) 




(12) 



In this case, g(L) remains constant in the critical point and its general behavior is shown 
in Fig. l,a. The relation (12) was never proved but its validity can be expected by analogy 
with the usual phase transition theory. 

The solution ip n {r±) of the Cauchy problem for the quasi-lD Schroedinger equation 
with the initial conditions on the left edge allows decomposition 



,72" 



+ 



(13) 



1 There are also pure imaginary (3 S , which are inessential for us. The detailed form of solution for 
Xm.m' (n) is given by Eq.44 in pjj and demonstrates equality of localization lengths in transversal and 
longitudional directions. 

2 Analogous calculations for E ^ give = 2\/4 — E 2 in accordance with pp. 

In comparison with [7] , it is based on more general motivation and is not restricted by one-parameter 
scaling hypothesis. 
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Figure 1: a — Typical behavior of g(L) in the case of one-parameter scaling: different 
curves for W > W c or W < W c can be made coinciding by a scale transformation, b - 
Behavior of g{L) according to Eq. 9: crude violation of scaling is evident. 
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where r± is the transverse coordinate (like m in (1)), the quantities hffl(r±) have no system- 
atic growth in n, while the Lyapunov exponents 7 S tend to the constant values in the large 
n limit. According to Mott [Sj, eigenf unctions of a quasi- ID system can be constructed by 
matching two solutions of the type (13) increasing from two opposite edges of the system. 
The tails of the eigenfunction will be determined by the minimal positive Lyapunov ex- 
ponent 'jmin and £i£> can be estimated as l/jmin] this estimation in combination with the 
scaling relation (12) is a basis for the most popular of contemporary numerical algorithms 
(see ffj and a review article [9]). Due to evident analogy between (5) and (13) we can 
refer (as suggested in [T]) to the exponents (3 S as generalized Lyapunov exponents. It is 
easy to show [3] , that positive (3 S and positive 7 S are in one to one correspondence: a term 
containing f3 s is an averaged square of the term containing 7 S . The usual correspondence 
between the mean square and the most probable value gives inequality (3 S > 2^ s [3]. 

Though no rigorous relation exists between (3 min and •y m i n , these quantities are very 
close from the physical viewpoint. Indeed, according to [3j: (a) inequality (3 m i n > 2^ min 
can be rigorously proven; (b) the order of magnitude relation /3 min ~ 7 min takes place in 
the typical physical situation^; (c) f3 min and 7 min are practically equivalent from viewpoint 
of one-parameter scaling philosophy, and relations 



can be postulated on the same level of rigorousness. 

If the correlation length £id is estimated as l//3 m j n , then the behavior of g(L) determined 
by Eq. 9 has a form presented in Fig. l,b; one can see an essential difference from a typical 
scaling situation (Fig. l,a). Since there is no growth of g(L) for large L, the state with 
long-range order (i.e. metallic phase) is absent, in accordance with [10] . Exponential 
localization takes place for W > W c , while the behavior specific for a critical point is 
realized in all range W < W c . The latter situation corresponds to localization with the 
divergent correlation length and probably should be interpreted as power law localization. 
The transition at W = W c is of the Kosterlitz-Thouless type and should not be mixed with 
the usual Anderson transition. This result can explain the observable 2D metal-insulator 
transition [TT]0 and does not imply a serious revision in the weak localization region. 

Violation of scaling relation (12) is clear in Fig. l,b since g(L) is not constant for W = 
W c . It is still more evident for W < W c , when different curves have different constant 
limits for L — > oo and surely cannot be made coinciding by a scale transformation. 

4 This relation follows from existence of the log-normal distribution and the fact that fluctuations of 
7 S are of the same order (or less) as their mean value. These properties can be proved in the limits of 
weak and strong disorder and are confirmed by extensive numerical studies in the intermediate region (see 
references in [51 fT^]). 

5 According to Last and Thouless [T5] , the hopping conductivity for power-localized states goes to zero 
for T — > as a power of T . Consequently, decrease of resistivity with growth of T is not so quick as 
for exponentially localized states and can be changed to increase by other effects (for example, by T 2 
contribution from electron-electron scattering). Such picture is close to observable in experiment [14]. 




and 




(14) 
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One may suspect that these results are related with our estimation of as 1/ j3 min , since 
the quantities (3 min and 7 mm can be essentially different. In fact, inequality (3 min > 2j min 
is sufficient for existence of phase transition. Indeed, one can see from this inequality that 
7mm — ► for L — > oo, if W < W c ] on the other hand, for large W existence of exponential 
localization is beyond any doubt and finiteness of 7 mm is evident; it can be rigorously 
proved in the large W limit [T3j. Of course, the upper bound for 7 min does not forbid it 
to decrease more rapidly than l/L, as it should be for a true metallic state. But such 
possibility is reliably excluded by numerical studies. 

The latter fact, in combination with inequality /3 mm > 27 mm , is sufficient to claim 
violation of scaling for 7 mm . Indeed, for W = Wo < W c the scaling parameter \/j min L 
cannot increase and is bounded from below by constant 2/f3 min L. If we take W\ < W 
such that PminiWi) < 2j min (W ) , then 7 mm (H / i) < 7mm(Wo) and the constant limit of the 
scaling parameter l/j min L cannot be the same for different W: we return to the picture 
presented in Fig.l,b. As a result, substitution of 7 mm for /3 mm does not lead to qualitative 
changes in the presented picture and inequality /3 mm > 2^ min is sufficient for the most 
responsible statements. 

One can see that interpretation in terms of finite size scaling leads to unambiguous 
conclusion on existence of the 2D phase transition. This statement does not contradict 
to numerical results, if the raw data are considered [3]. The opposite conclusion made by 
numerical researchers is based on interpretation in terms of one-parameter scaling, which 
is inadmissible here. Recent numerical results clearly demonstrate (see Fig. 37 in [9]) that 
possibility of the 2D phase transition cannot be rejected on numerical grounds. 

However, we do not consider the conventional variant of finite size scaling as indis- 
putable. Calculation of Lyapunov exponents is not equivalent to diagonalization of the 
Hamiltonian: for example, statistical independence of i\) n ,m an d V n ,m' is valid in the first 
but not the second case, and one can suspect oversimplification of the problem. For reli- 
able estimation of £i£> one needs detailed study of the coefficients A4 in Eq.13 appearing in 
Mott's construction for eigenf unctions. In the general case, relation £ 1D ~ l/7 mm can be 
violated and some effective exponent j e ff should be used instead 7 mm [3]. If the scaling 
relation of type (14) is postulated for Jeff, then its dependence on parameters is determined 
by scaling itself. Analysis shows possibility of two variants: (a) the 2D phase transition is 
removed; (b) the 2D transition remains, though behavior of the correlation length becomes 
different. One can suggest that both possibilities are realized in different models. 

The described approach can be generalized to higher dimensionality [TBI C2] • Unfortu- 
nately, the analytical results for the critical disorder appear to be in essential contradiction 
with corresponding numerical results. The interpretations of this fact can be different, but 
in any case it is related with crude violation of one-parameter scaling [T2] . 

3. Approach by Kuzovkov et al JTJ 

Let us discuss the difference of approach presented in Sec. 2 from one suggested by 
Kuzovkov et al [JJ, US]. The initial system of equations (4) and its higher dimensional 
analogue (Eq. 5 in [12]) coincide with those used in [HUE]- However, the quantity z mm i{n) 
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was not introduced in [HUB] and its role was played by y m 'm( n )- As a result, the system of 
equations had no complete difference form and could not be solved in the natural manner 
with evaluation of full spectrum of exponents (3 S . Instead, the authors of [HE] used the Z- 
transform and introduced the so called filter function H(z). The latter, as was demonstrated 
on simple examples, has the poles corresponding to eigenvalues of the transfer matrix, 
and, in principle, this approach allows to find the full spectrum of generalized Lyapunov 
exponents. However, some problems arise in the practical realization of this scheme: 

(a) The solution was obtained only in the thermodynamic limit L —>■ oo, and the 
problems exist with interpretation of results (see below). 

(b) Averaging over translations in the transversal direction was used for simplification 
of the problem. This procedure is disputable since it can lead to elimination of some singu- 
larities of H(z). Indeed, translational invariance for the solution (7) of Eq. 4 and its higher 
dimensional analog (Eq. 10 in [12]) is absent and averaging over translations eliminates all 
terms with transverse momentum p ^ 0. Analogously, averaging of the squared solution 
eliminates all terms with p ^ 0, G/2 where G is a vector of a reciprocal lattice corre- 
sponding to the main diagonal of the Brillouin zone. We can suggest by comparing results 
(though cannot prove it rigorously) that the filter function H(z) is analogous to i\)\ m and 
only terms with p = and p = G/2 are essential for it. Fortunately for the authors of 
[U [16], a condition p = G/2 corresponds to the minimal exponent f3 m i n for d — 2 and 
d > 4, so the critical values a' (corresponding to W c in P, [12]) were found correctly for 
these cases. In the case d = 3, the exponent f3 m i n does not correspond to p = G/2 and 
a zero value for the critical disorder was obtained in [TU] for the band center E = 0, in a 
striking contrast with [T2] . 

Attempt to justify this point made in the paper [17] is not convincing, in our opinion. 
In fact, the results reduce to the statement that averaging procedure does not eliminate 
the maximum exponent /3 ma x, corresponding to p = 0. This statement is surely correct 
but it has no relation to physics of the problem. 

(c) In the general case, function H(z) can have not only poles corresponding to eigen- 
values of the transfer matrix but also another singularities, which are physically irrelevant. 
In our opinion, it is an origin of the second critical point <jq obtained in [JJ [16] for higher 
dimensionality; we see no evidence for it in the spectrum of f3 s . Correspondingly, we see no 
evidence of a special role of dimensionality d = 6, which is surely absent in the exact field 
theory approach [IS] . 

The papers [TJ [16] were formulated in the engineer's language (using the concepts of 
signals, filters etc.), which has no direct relation to the Anderson transition. The authors 
relate the extended states with a stable filter and the localized states with unstable filter, 
i.e. introduce their own localization criterion whose correspondence with conventional one 
was never studied. 

The limit L — > oo was taken in [TJ [16] from the very beginning, and the finite-size 
scaling approach could not be used for interpretation of results. The authors interpreted 
the Anderson transition as being of the first order, considering two branches of the filter 
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function H(z) as two different phases existing simultaneously; such interpretation does not 
have clear physical sense. 

Possibility of the first order transition is in conflict with the old ideas by Mott [8]: 
according to the Poincare theorem, a small change in the energy or the disorder strength 
induces small changes of the wave functions, and hence the state of the system changes 
continuously^. It is essential, that the Poincare theorem is valid only for a finite system 
and allows existence of the mobility edge, if the localization radius has divergency in it. 
Existence of the divergent length is not characteristic for the first order transition. 

The possibility of power- law localization was also mentioned in [U [16], but in the con- 
text, which is completely different from that in [31 [T2]. 



4. Approach by Markos et al |2J 

Markos et al [2] do not follow the natural procedure of Sec. 2, but advance more com- 
plicated approach. They rewrite Eq. 1 using the transfer matrix and construct the direct 
product of two such matrices. After averaging, they arrive to a linear system of equations 
determined by a matrix 



T 



( W 2 + D ® D -A)®1 -1®A) 1(8)1 ^ 

D (8)1 -1®1 

1®A) -1(8)1 

V i(8)i o o o J 



(15) 



where D = E — H Q and H Q is the Hamiltonian of the n-th slice for a pure system, 1 is 
the unit matrix of the size L x L. In principle, this approach is equivalent to that of Sec. 2 
but appears to be practically untractable due to sophisticated matrix constructions. The 
eigenvalues A = exp(ig) of the matrix (15) are declared to be determined by equation 

2cos2g-2/t i /t j cosg+(^ + K J 2 -2) = 2W 2 isinq (16) 

where Ki = E — 2 cos Pi and p^ are allowed values of the transverse momentum p. Equation 
(16) surely not coincides with the corresponding Eq. 9 in Sec. 2. The main difference is the 
absence of functions with the argument qL, which inevitably arise due to the boundary 
conditions and can be absent only in trivial cases. One can suspect, that this difference 
is related with improper treatment of the disorder term (like W 2 1 x 1 in (15)), without 
which the problem is trivial. This term is local in m — wl (see (4)) and not diagonal in the 
momentum representation, which is seemingly used in (16). It looks likely that the local 
nature of this term was neglected and it was replaced by a suitable constant. We can try 
such thing for the system (4), replacing 5 mjm > by unity. Then (4) is solved trivially 

rr ,( rl \ — ^pipm+ip'm' J3n ( \ _ ipm+ip'm' fin I \ _ ipm+ip'm' fin 

J 'm,,m'V"7 J/C c ) yrn,m' \' J i/ c c j '^m,m'\ 11 '} ^ c c j 

(17) 



6 It does not exclude that certain quantities can display a jump-like behavior; the conductivity, according 
to Mott, belongs to such quantities. 
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where allowed values for p and p' (2ns /L, s = 0, 1, . . . , L — 1) are determined by the 
boundary conditions and the quantities x, y, z satisfy the equation 



W 2 + 4 cos p cos p' 
2 cosp 
2 cosp' 



2 sinh p 



2 cosp 


2 cosp' \ 










1 


i J 





(1« 



The determinant vanishes under condition 



2 cosh 2/3 - 2/c« / cosh (3 + (k 2 + k' 2 - 2) = 2W 2 sinh f3 



(19) 



where k = —2 cosp, k' = —2 cosp'. In the case E = (when Eq. 4 holds) this equation 
is identical to (16), if correspondence (3 = iq is taken into account. We see that, indeed, 
the disorder term in (15) was treated inadequately and its local nature was neglected. 
Physically, equation (16) corresponds not to the true Anderson model but to its degenerate 
version when site energies V n ^ m are independent of m. 

In fact, the error is present already in the matrix (15). To obtain the term with disorder, 
one needs to produce averaging of the kind 











( VxVx 








\ 








°) = 





V1V2 












9(o 


v 2 J 








v 2 v l 






V 







w 2 



V 2 V 2 ) 



( 1 




V 













1 / 



(20) 

and the result cannot be represented as W 2 1 ® 1 . The latter form is valid for V\ = V 2 in 
Eq. 20, while in the general case it corresponds to the model (1) with V n ^ m being independent 
of m, in accordance with the previous analysis. Such model is of no interest and all 
conclusions made in [2j are irrelevant for the problem under consideration. 



5. Conclusion 

The following conclusions can be made from our analysis of generalized Lyapunov ex- 
ponents: 

(a) A spectrum of f3 s allows complete analytical investigation. 

(b) Interpretation of results in terms of conventional variant of finite-size scaling leads 
to unambiguous conclusion on existence of the 2D transition of the Kosterlitz-Thouless 
type and crude violation of one-parameter scaling in any dimensionality. 

(c) Attempt to restore one-parameter scaling [TU] by replacement 7 min by j e ff eliminates 
the 2D transition from roughly half of models. 

(d) Results for the critical disorder obtained by Kuzovkov et al [U [16] are correct for 
d = 2 and d > 4 (not for d = 3), but their interpretation is not satisfactory. 

(e) Qualitatively different conclusions by Markos et al [2] are based on improper calcu- 
lation of the matrix product. 
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(f ) Above conclusions do not contradict numerical results, if the raw data are considered; 
interpretation of the latter in terms of one-parameter scaling is inadmissible. 

This work is partially supported by RFBR (grant 06-02-17541). 
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